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Abstract 

We study massless QEDi with N flavors using path integrals. We 
identify the sector that is generated by the N 2 classically conserved 
vector currents. One of them (the U(l) current) creates a massive par- 
ticle, while the others create massless ones. We show that the mass 
spectrum obeys a Witten-Veneziano type formula. Two theorems on 
n-point functions clarify the structure of the Hilbert space. Evaluation 
of the Fredenhagen-Marcu order parameter indicates that a confining 
force exists only between charges that are integer multiples of ±Ne, 
whereas charges that are nonzero mod(iV) screen their confining forces 
and lead to non-vacuum sectors. Finally we identify operators that 
violate clustering, and decompose the theory into clustering 8 vacua. 



* Also at Institut f. Theoretische Physik d. Universitat Graz, Austria 



1 Introduction 



The Schwinger model is a laboratory to study many important aspects of more 
realistic models of particle physics such as QCD. It shows confinement, mass 
generation of its would-be Goldstone particle via the axial anomaly - thereby 
solving its U(l) problem, it also has topological sectors and a vacuum angle 
9 dual to them. While these features are well known and understood in the 
case of Schwinger's original one-flavor model, in a model with several flavors 
the situtation is much less trivial and has not been investigated in much detail 
so far. This iV-flavor Schwinger model has in fact more features in common 
with four dimensional QCD: On the classical level it has a symmetry group 
U(N) L x U(N) R that is broken down by the anomaly to SU(N) L x SU(N) R x 
U(l) v just like in QCD. 

In this paper we study this model using the euclidean functional integral. 
The boundary conditions we use (vector potentials decaying at infinity suffi- 
ciently fast to be square integrable) correspond to zero topological charge for 
the gauge field. This leads to an infinite volume vacuum state enjoying chiral 
invariance, but violating the cluster decomposition property. It turns out that 
by explicit construction this vacuum state can be decomposed into a mixture 
of clustering states labeled by an angle 9, the so-called vacuum angle, which 
is conventionally introduced as a dual variable to the topological charge (first 
Chern number) of the gauge field. 

We compute the correlation functions of the vector /pseudovector currents 
as well as the scalar /pseudoscalar densities and find that the theory describes 
a multiplet of one massive and N 2 — 1 massless pseudoscalar 'mesons'. There 
is a relation between the masses that looks exactly like the so-called Witten- 
Veneziano formula for the pseudoscalar meson masses in QCD; in fact it allows 
to give a precise meaning to this formula with proper attention to the short 
distance fluctuations of the topological density that play an essential role. 

The model also reveals an interesting superselection structure. Computa- 
tion of the so-called Fredenhagen-Marcu confinement parameter suggests that 
there are superselection sectors labeled by a charge that is defined as an integer 
mod(iV). This means that the confining Coulomb potential does not prevent 
operators of fermion number ^ mod(iV) from creating new charged sectors; 
their electric charge seems to be 'screened'. Of course QCD is not expected to 
share this property. 

The model is still trivial in the sense that there is no real interaction be- 
tween its particles, but it is not Gaussian like the one-flavor model. Its vac- 
uum sector turns out to be a tensor product of the Fock space of a free massive 
scalar particle with the representation space of the SU (N) L x SU (N) R current 
(Kac-Moody) algebra induced by N massless Dirac fermions. 

We would like to make a remark about the level of rigor of this paper: Our 
aim is to elucidate the physical content by explicit calculation. Therefore we 
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do not spend time to justify all the manipulations involved in deriving the 
results with full mathematical rigor, even though it would not be too difficult 
to do so. In |IJ it is sketched how this is to be done for the one-flavor model; 
the N— flavor models does not pose any fundamentally new problems in this 
respect. 

There are some open problems left: For instance it should be expected 
that the subspace created form the vacuum by the scalar /pseudovector cur- 
rents is identical to the one created by the vector /pseudoscalar densities; we 
did not attempt to prove that. Likewise the structure of superselection sec- 
tors suggested by the Fredenhagen-Marcu parameter has not been established 
rigorously, because that parameter by its nature has only suggestive value. 

The paper is organized as follows: In the next section we define the model 
and introduce our notations. Furthermore we derive some basic expressions 
for generating functionals and Green's functions. In Section 3 we investigate 
the sector that is generated by vector currents. The mass spectrum that we 
obtain fits a Witten-Veneziano type formula, which we discuss in Section 4. To 
get more insight into the subspace generated by the vector currents we derive 
general expressions for their n-point functions in section 5. The next section 
is dedicated to the problem of confinement. We evaluate the Fredenhagen- 
Marcu order parameter for a e + e~ system and higher products of fermion 
fields. Finally we decompose the model into clustering theories and thereby 
obtain the vacuum angle. 



2 Setup 

We use the following representation of the 2-d Euclidean 7-matrices. 

7l= (? l) ' 72= (. ~o) ' 75 = 27271 = (J _°i) (!) 

They obey the commutation relations {j^ju} = 25^ ; v — 1,2,5. The 
Green's function G(x,y;A) for fermions in a background field A obeys the 
equation 

luidp — ieA^)G(x, y;A) — 6(x,y) . (2) 

In two dimensions it can be solved explicitely ||, provided the background 
field satisfies some mild regularity and falloff conditions; the solution is 

G(x, y; A) = e ^)-^y)] G ^ x _ y ) (3) 



where 



<&(x) = - I d 2 zD(x - z 



d^A^z) + i^Efu, (z) . (4) 



G (x) is the Green's function for zero background field 



G\x) = f 7 ^ . (5) 
2tt x z 
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D(x) denotes the Green's function for free massless bosons (—AD(x) = S(x)). 
The action for the iV-flavor Schwinger model reads 



S[il>,1>,A] = S g [A]+S f [ii>,1>,A] , (6) 

with 

1 



S g [A] = J d 2 x 



-F^(x)F^(x) + -\(d,A,(x)Y 



^QA). (7) 



The kernel Q of the quadratic form we introduced is given by 

Q = (1 — A)<9®<9 — A . We work in Landau gauge, which means A — ► oo after 

the integration over the gauge field. 

The fermionic part of the action is a generalization of the one flavor model 



- N r ■ 

S f ty,1>,A] = £ / d 2 x^ a \x)^[d,-teA,(x)}^ 

a=l J 

ip^ a \ ip^ are 2-component spinors with a flavor index a = 1, 2, ...N. It is useful 
to introduce sources for the fermions and to define the generating functional 



(9) 

Vacuum expectation values can be obtained as functional derivatives with 
respect to r],r]. Performing formally the Berezin integral over the Grassmann 
variables leads to 

Z[r),rj] = JidA]e- s ^detl0-teA] N e^J d2xd2y ^ . (10) 

Note that the formal expression det[$—ie is only defined when an ultraviolet 
and infrared cutoff (for instance a finite space-time lattice) is introduced. The 
determinant can then be normalized to 1 for e = 0, by replacing it with 
det[l — eK(A)] where K(A) = i J]L<f)~ . In two dimensions this determinant can 
be computed explicitely 0, using the idea of regularized fermion determinants 
(for a review see e.g. 0, for the QED 2 case [[I]). If we assume that the vector 
potential A satisfies some mild regularity and falloff conditions at infinity to 
make it square integrable, the answer is 

det[l - ei^(v4)] = e-i^il^Hi , (11) 
where A^ is the transverse part of the gauge field 



Al = A,-^A u . (12) 
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We combine the logarithm of the determinant and the gauge field action in 
one common quadratic form for A^. This gives rise to a Gaussian measure 
dfi c [A] = \[dA\ exp(— \{A, C~ l A)) for the gauge field with covariance C de- 
fined by 

(A, C~ l A) := (A, QA) + e 2 - || A T \\\ . (13) 

7T 

A few lines of algebra lead to 

C-l = ( - A + e 2 ^) (V - dj §) ~ d,8A . (14) 
Inverting and taking the limit A — > oo gives 

7T 

Besides this form of the measure for the gauge field, we compute the measure 
for a field (p which is a linear functional of the gauge field. This leads to a very 
simple dependence of the propagator on the gauge field. G(x, y; A) depends on 
$L4] via the exponential exp(ie[$(x) — $(y)]) (see equation (3)). <£> in general 
depends on the gauge field through the terms d^A^ and e ilv d il A v ^ but the 
first term vanishes (with probability 1) in Landau gauge so that $ becomes 
a function of e^ u d^A v alone. To get the spinor structure of the propagator 
explicitly, we define (after dropping the d^A^ term in $ ) 

$(x) =: ij 5 (p(x) , (16) 

where 

<p(x) = - J d 2 zD(x - z)e ilv d iX A v {z) . (17) 

This implies 

V? = £j ^A u . (18) 

Using this equation one obtains the Gaussian measure d^g[ip\ for the field ip, 
induced by the measure for the gauge field with covariance C given by 

C = . (19) 

(-A + 7Vf)A 

Expressed in terms of ip the propagator reads 

G(x,y; <p) = e~ e ^ {x) -^ y)] G°{x - y) = 



1 1 / e -eMx)-v(y)] ( 5 _ 

2tt (a: - y) 2 \ e +e ^~^ (x - y) 
where we introduced the complex coordinate 



(20) 



x = x\ + ix2 ■ (21) 
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3 Vector currents 



In this section we derive the 2-point functions of the vector/pseudovector cur- 
rents and obtain thereby information about the particle spectrum of the theory. 
Note that the axial vector currents are simply linear combinations of the vector 
currents due to the fact that in two dimensions 757^ = is^v 



3.1 Bosonizing a Cartan subalgebra of U(N) 

It was noticed in ||] that the U(l) current together with the other elements 
of a Cartan subalgebra of U(N)fi avor can be bosonized in terms of free scalar 
fields. This result is the first step in our analysis of the particle structure of 
the model and therefore we reproduce this result in terms of the functional 
integral formalism. 

Define vector currents 

j®(x) :=^ {b \x)^ b \x) ,b = 1,2,. ..N. (22) 

In principle one would have to regularize these currents by e.g. a point splitting 
prescription, but we can use the fact that the fermion current couples to an 
external field a M the same way as to the quantized field and thereby avoid all 
problems due to short distance singularities. We couple the currents to source 
fields a^\x) obeying 9^' = 0, and evaluate the characteristic functional 
using the explicit form of the determinant 

A? 



6=1 



-(A,A) = e -#Ej a(i,) Hi e |(A,CA) ^ 



e 2 N 



.fELctf^)) = y [ d A}e- s ^ J] det[l - eK(A + a «)" 

e -#Ejl« w lll Jdixc[A\e 

A, = -E< } - (23) 
n 6=1 

After some reordering of the terms in the last equation one ends up with the 
following expression for the characteristic functional 

where M is given by (using matrix notation in flavor space) 

-e 2 /vr 1 



-A + e 2 ^)A -A + e 2 ^' 
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R 



(N-l -1 
-1 N-l 
-1 



-1 \ 



-1 



(25) 



\ -1 . —1 N — 1 J 

The numerical matrix R can be diagonalized by an orthogonal matrix U given 

by U T = (r^\r^ 2 \ r^) where the are the normalized eigenvectors 

of R. The corresponding eigenvalues are 0,N,N, N, implying URU T = 

diag(0, N, ....N). This allows to diagonalize the covariance M. Define 



K := UMU 1 



-A+e 



2 N 



1 

-A 



V 



1 

— ) 



(26) 



We now choose the sources in flavor space proportional to one of the 
eigenvectors r®. This corresponds to a change of the basis in flavor space and 
allows to express expectation values of the new vector currents in terms of free 
bosons. Define 

JP := ^if Cty ? (27) 

where the N x N matrices H"> are generators of a Cartan subalgebra of 
U(N) 

flavor 



N — 1 + (N — 1) 
1 



N 



diag(l, 1, ....1, —N + 1) , 



^N-2 + (N- 2) 2 
1 



diag(l,l, 1, -TV + 2,0), 



#W = -J= diag(l, -1,0,0, 0). 



(28) 



The diagonals of these generators are the eigenvectors of M. Using equation 
(24) one easily finds the generating functional for the 

^-^[e^d^ ,K U £ p<T d p b^ _ ^ie-^^ v d vV ><- l \b ( ^\^ 

where bty is a source for the currents JP. The index I now labels the 'meson' 
fields (27). (•) B denotes the Gaussian expectation values of an iV-tuple of scalar 
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fields (p( l \l = 1,2,...N with covariance K given by equation (26). Hence we 
can identify 

JP = 4^<V Z) • (30) 
The last equation is the bosonization formula of a Cartan subalgebra of U(N) 



3.2 More meson fields 

In analogy to the construction of meson states in QCD, one can define vector 
currents for all the generators of U(N). A convenient basis of the Lie algebra 
of U(N) is given by the N(N — l)/2 generators of the form 



\ 



1 

V2 



(31) 



and N(N - l)/2 of the form 



1 

V2 



(32) 



V / 

The corresponding vector currents are 

Jjp(x) :=^(x)H il) tfj(x). (33) 

Again no point splitting has to be taken into account, since only different 
flavors (which cannot contract) sit at one space-time point and therefore there 
is no short distance singularity. 

First we notice that all N 2 vector currents generate orthogonal states 

(jP(x)jP(y)) = 5 n , y) Z = l,2,...iV 2 , (34) 

where JF is the two point function. For the Cartan currents (34) follows 
dircectly from the bosonization. To prove it for the set of all iV 2 currents 
one has to take functional derivatives of the generating functional (10) with 
respect to the fermion sources 77, 77. If I 7^ I', either the different flavors do not 
contract entirely, or terms with opposite sign cancel. 
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In the case I — I', I — N + 1, N + 2, ...N 2 , functional derivation leads to 
(take e.g. \i — v — 1) 

Fix = ~ J dfi <5 [ip]{G 2 i{x, y; ip)G 21 (y, x; if) + G 12 {x, y; (f)G 12 (y, x; </?)} . (35) 

Using the explicit form (20) for G, one immediately sees that the exponentials 
involving <p cancel. Integration over (p simply gives a factor 1. The same is 
true for arbitrary fi, v. Hence the two point function is the same as for free, 
massless fermions. It can be expressed in terms of derivatives of the propagator 
of a free massless boson, giving rise to the same expression as was obtained 
for the Cartan subalgebra, i.e. for I = 2,3,...N. Putting things together we 
conclude for the two-point functions of the vector currents 

(j«\x)jP(y)) = 5 IV h w dx p e va dy a {<p®{x)<p®{y)) B , (36) 

where the scalar field <p^ has mass m = and (pV\l = 2,3,...N 2 are 

massless. In addition we obtained explicit expressions for arbitrary n-point 
functions of the Cartan currents (/ = 1,2, ...N) via the bosonization prescrip- 
tion (30). In Section 5 we will obtain the n-point functions for all the currents. 

4 Mass spectrum and Witten-Veneziano type 
formulas 

In 1979 Witten || proposed a formula that relates the mass of the r/ meson 
to the topological susceptibility of quarkless QCD. The purpose of this was to 
show how QCD solves its U(l) problem, without having to bring instantons 
into the discussion. The derivation given by Witten has some problems (terms 
of equal sign are cancelled against each other) as was pointed out by Seiler and 
Stamatescu 0. But it was also noted there that for the one flavor Schwinger 
model the formula is nevertheless true, if the topological susceptibility is inter- 
preted apppropriately, with proper attention to contact terms (the fact that 
therefore the correct treatment of short distance fluctuations of the topologi- 
cal charge density is essential makes of course the application to lattice gauge 
theories cumbersome, to say the least). 

The N flavor Schwinger model is a convenient laboratory to study the U (1) 
problem in a more general setting, much more reminiscent of QCD, because 
it has both flavor changing axial currents that are conserved and a flavor non- 
changing one whose conservation is violated by the axial anomaly. The mass 
spectrum we obtained in the last section contains both massless Goldstone 
particles related to the conserved axial currents and a massive particle related 
to the anomalous U(l) axial current. So it has the features expected to be 



present in the meson sector of QCD in the chiral limit. Hence it is natural to 
check if our model satisfies a formula of the Witten-Veneziano type. 

Here we cite the more general formulation of Veneziano which re- 

duces to Witten's form when setting all meson masses (except m v ') equal zero. 
It is given by 

m v> ~ 2 m v ~ o m -° = ~~n x ■ ( 37 ) 

J 7V 

f n denotes the meson decay constant. % has a more subtle meaning. The 
interpretation as topological susceptibility faces some problems, as was pointed 
out by Seiler and Stamatescu ||. But a formula of this type can be derived if 
X is interpreted as the contact term that occurs in the spectral representation 
of the two point function of the topological charge density (compare also the 
appendix in ||, and the derivation given by Smit and Vink ||). Here we use 
the form that was obtained by Seiler and Stamatescu for an arbitrary number 
of flavors N and vanishing quark masses, which implies = m^o = 0. It 
reads 

P(0) denotes the announced contact term. (38) contains the rj' decay constant 
frf, which shows up also in Witten's derivation, and is replaced later by to 
which it is approximately equal. 

In QED 2 the topological charge density q{x) is given by 

q(x) = ^F l2 (x) . (39) 

The topological susceptibility is defined to be 

e 2 r e 2 



Xto P = J (q(x)q(0))dx = J (F 12 (x)F 12 (0)}dx = ^GW(O) , (40) 

where Gff denotes the F 12 propagator, and Gff(Q) h's Fourier transform at 
zero momentum. Since 

F 12 (x) = e^A u (x) , (41) 
the propagator Gff is given by 

Gff — —^^luduCppEpada , (42) 

where C is the gauge field propagator (15). Inserting C and transforming to 
momentum space we obtain 

g„(p)=-^W=i- -J^=i-r*« . (43) 

V + e P + e J° P + A 4 
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In the last step we made the spectral integral explicit. One nicely sees that 
the spectral measure 

dp(/i 2 ) = 5(/i 2 - m 2 )dp 2 (44) 

is 'dominated' by the contribution of the rj', which should of course be identified 
with our massive particle. Hence we have 

m i = m = e — . (45) 

' 71 

From (43) one immediately reads off the contact term -P(O) in the spectral 
decomposition of Xto P 

P(0) = ^ • (46) 
The last missing ingredient is the decay constant /,,/. It is defined by 

f lf =m- 2 (0\d lt jV>\T,'). (47) 
The anomaly equation for the current jjfi = ip^'y^ip^ reads (see eg. |jj ) 

dtfP = 2? • (48) 

Using the definition of the U(l) current (27), (28) one obtains the anomaly 
equation^] 

d^J {1) = , (49) 

which we insert in (47) to end up with 

/„/ = m~ 2 2\fN(0 \q\rf) = m~ 2 2\^N(0 \ q —F 12 | 0) . (50) 

m 

In the last step we generated the \t]') state as Z^Fi 2 \0) with the normalization 
condition 

Gwn' = ~t o + contact term , (51) 

■p 1 + m l 

giving rise to Z = — We end up with 

U = i\ ■ (52) 



Insertion in (38) of the quantities obtained gives an identity. This explicit 
computation shows that eq. (38) holds in massless QED 2 with N flavors. 
But one can also interprete the result as a verification of the original form 
of the Witten-Veneziano formula, because the topological susceptibility of the 
quenched theory reduces to the contact term. It is not true, however, that the 
topological susceptibility appearing in the formula expresses only a property 
of the long distance fluctuations of the topological density. 

1 Note that in our definition of the 17(1) current there is an extra factor 1/vJV which 
modifies (47), (49) by this factor, compared to the usual notation. 
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5 N-point functions 



To gain more insight into the sector generated by the vector currents we com- 
pute explicit expressions for the connected n-point functions C n . Consider- 
ing fully connected correlations has the advantage of avoiding contractions of 
fermions at the same point. Therefore we can avoid using test functions as 
would be required in principle by the distributional nature of the fields. Let 

Cn+i == (j^mjSM) J S: } M) c = 

^a bo^aibi ^a n b n (7/uo)aoA)(7Mi)ai/3i {l/J.n)a n f3 n 
o.i,bi oti,Pi 

(^(^^(xo^^xo^^xo ^::\x n )4:\ Xn )) c . ^ 

Nonvanishing contributions occur only if the color indices can form a closed 
chain (e.g. bo — ai, b± — a 2 , ■ ■■b n -i = a n , b n = ao). The corresponding factor is 
simply the trace over the flavor matrices H®. To find all possible contributions 
one has to sum over all permutations it, keeping the first term fixed. 

C n+l = - £ Tr[HM H M HM] E(7 w ,)«oA,(7w)«ift (7,J^„ 

7r(l,2,...n) ai,/3i 



I' 



^c[#A« l( i ) ( 3; o,^(i)^)G ft(1)ai(2) (^(i),^(2)i <p) Gp w(n)ao (x n{n) ,x ; (p) 

(54) 

Since in the chosen representation the 7 matrices and the propagator G(x, y; (p) 
have only off-diagonal entries (cf. equation(20)), we find the following chain 
of implications for e.g. fa — 1 

fa = 1 =>- ^(1) = 2 =>- ^(1) = 1 =>- a %{ 2) = 2 f3 n ( n ) = 1 =>• a = 2. (55) 

When starting with fa = 2 one ends up with the opposite result of all fa = 2 
and all ccj = 1. Besides those two no other nonvanishing terms contribute. 
Again one finds that all dependence on ip cancels, and only free propagators 
G° remain. Using G°(x) = -^p^r- and making use of the complex notation 
(21) we obtain 

C n+1 = - 1 - Y Tr\ H^H^ 

(27r) ri+1 ^ L 



7r(l,2...n) 



( n 11 1 1 

{lI(7 W )2i- — + cc \ , (56) 

where c.c. denotes complex conjugation. 

Using this basic formula, we construct for N > 2 and arbitrary n fully 
connected n-point functions that do not vanish; in other words we show that 
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the currents are not Gaussian. For simplicity we give the construction for the 
case of N = 2. Since for arbitray N > 2 generators with the same commutation 
relations exist, it is obvious how to generalize the construction to general N. 
The generators we are using are the Pauli matrices up to a normalization factor. 
To distinguish them from arbitrary generators we denote the special set 
of matrices needed in the construction 



(i) - 



1 



T 



Define 



(° 1 



-(2) 



1 

71 



-i 

1 



-(3) 



1 



(i o 



V2\0 



F"n+i(xQ, Xi, ...x n ) : — 

1 1 

Xq — 5tt(1) ^tt(I) — ^tt(2) 

where the special set of r matrices is given by 



E Tr 

7r(l,2,...n) 



t ( 1 0) t (Itt(1)) 



1 



Xftlji) Xq 



(57) 



(5f 



{r 



(io) T (h) 



-(In) 



f ^(2)^(3)^(3)^ ^(3)^(1)} {orn + i =2m + l (n + 1 odd) 
| { r (i) , r (3) ? r (3) ) _ _ _ r (3) ? r (i) } for n + i = 2m + 2 (n + 1 even) 

Theorem 1 : 

For arbitrary m > 1 : 

-P2m+1 ^0 , I*2m+2 ^0 . 

Proof: We prove the statement by induction. 

1) F 3 ^0 and F 4 ^0 can be checked easily. 

2) Assume F 2m +i ^0 and F 2m +2 ^0 . 

3) We have to show F 2m+ 3 ^0 and -F 2m+ 4 ^0. First we check -F 2m+3 . 
we abbreviate 2m + 3 := n + 1. The trick is to consider F n+1 (xo, • 
a function of x , and to compute the corresponding residues 0. 



(59) 



(60) 



Again 
■Xn. ) as 



Res £o [F n+ i,x 



E Tr 

7r(l,2,...n),7r(l) = l 



r (2) r (3) r (7r(2))_ _ Wn)) 



1 



1 



1 



- E Tr 

7r'(l,2,...n),7r'(n)=l 



r (2) r ( 7 r'(l))___ r (7r'(r l -l)) r (3) 



x \ — 3+(2) a?7r(2) — ^7r(3) 2+(n) — ^1 
1 1 1 



Xi — av(i) 2V(i) — av( 2 ) av(n-i) — x\ 

„ (61) 

For a given permutation it choose the permutation n' such that 7r'(l) = 
7r(2),7r / (2) = 7r(3), ...ir'(n — 1) = vr(n), 7r'(n) = 7r(l) = 1. This gives for the 
residue 

1 1 1 



E 

(l,2,...n),7r(l)=l Xl X v(2) 3+(2) 3+(3) x ir(n) #1 



2 We acknowledge a useful discussion with Peter Breitenlohner on this point. 
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{Tr 



r (2) r (3) r M2)) r M3))___ r Mn)) 



Tr 



T (2) T ( 7 r(2)) T ( 7 r(3))__ T (7r(r l )) T (3) 



} • (62) 



Using the cycliticity of the trace and {t^ 2 \ t^} = 0, one finds that the second 
trace is the negative of the first one. Furthermore t^ 2 't^> = -^t^ . Looking 
at the definition F n for even n we conclude 



Res £o [F n+ i,x = x x ) = iV2F n \ 



.x r 



2771+2^1; •••2 ; 2m+2, 



(63) 



Since by assumption F 2m+2 ^0, Res[F 2m+3 ] ^ and hence -F 2m+3 does not 
vanish. The same trick can be applied to prove that this implies -F 2m +4 ^0. □ 

A n + 1-point function with the flavor content given by (59) is simply the 
real part of a multiple of F n+ i (compare(56)). Hence there exist nonvanishing, 
fully connected n-point functions for arbitrary n. 

This result implies that it is not possible to bosonize the whole set of all N 2 
vector currents J^> in terms of free bosons (this was the reason why Witten || 
introduced his nonabelian bosonization). This can be done only for a Cartan 
subalgebra where all generators commute, and the two traces in equation (62) 
cancel. The last observation allows to prove a second theorem. 
Theorem 2 : 

Any fully connected n-point function 



jSMj^( Xi )j^\x 2 ) j« : \x t 



(64) 



containing the U(l) current vanishes, with the exception of the two point func- 
tion (J^(x)jj) 1 \y)) c (without loss of generality we shifted the U(l) current in 
equation (64) to the first position). 

Proof : To prove the statement we use the same trick as for theorem 1. Again 
we consider the residues of functions F n+ i defined by 



Tr H^H^ 1 ^ 

7r(l,2,...n) 



Fn+l(Xo, Xi, ...X n ) :- 
1 



1 



1 



^0 ~~ x n{l) x tt(1) ~ ^7r(2) ^Tr(n) ~~ ^0 



(65) 



For the residue at x = x\ we obtain an equation equivalent to (61). After 
making the same choice for tt' one finds 



Res So [F n+1 ,x = xx] 



1 



7r(l,2,...n),7r(l) = l Xl X ^{ 2 ) X ^i 2 ) 



1 



— X 



t(3) 



X 



7r(n) X\ 



|Tr H (' 1 )ij( / T(2))...ff( / 7r(i)) — Tr i/( 1 )ij(^(2))...7f(^(n))_f/'( i l) | 



(66) 



Since H^ 1 ' commutes with all generators H^ l \ the two traces are the same and 

cancel. Using the same argument one can show that all residues at x±, x 2 , x n 

vanish. So F n+ i is analytic and bounded in the entire £q plane. By Liouville's 
theorem F n+ i is a constant, and the limit x — ► oo shows that this constant 
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is zero. Since the n + 1-point function defined in (64) is proportional to F n+ i, 
it has to vanish. □ 

Theorem 2 allows to prove the following proposition about the structure of 
the Hilbert space. 
Proposition : 

The Hilbert space 7i generated by the vector currents from the vacuum is the 
tensor product 

H = Hu(l) <8> Hmass=0 ■ (67) 

To prove this we need the connection between untruncated and fully connected 
n-point functions (see e.g. [|Kj). 

(01 0n)= E II<&1 ( 68 ) 

TTEVn PS7T 

where V n is the set of all partitions of {1,2, ...n}, n = {pi,P2, ■ ■■P\n\} denotes 
an element of V n , and {ii,i2, ■ ■■'i\p\} is an element p of it. An arbitrary n + k- 
point function (without loss of generality we write the U(l) currents first; 
li 7^ 1,2 = 1, 2, ...k) factorizes due to Theorem 2 



e n<^^)^ 

irGV n+k P^tt 



. X I 



"31 



e n<jgwjS' 



X ; 



TteV n PGvr 



>l 



e n (4^j4fw 



[JftMJgM Jff W) (^(i/O^CaO 4[ k \y k )) . (69) 

From this the tensor product structure of the Hilbert space follows easily. □ 

We observed in equations (54)- (56) that the dependence of the fully con- 
nected correlations on the gauge field cancels entirely. Hence the vector cur- 
rents in the 'massless' sector of the Hilbert space obey the same algebra as in a 
system of uncoupled fermions; this is the well-known level 1 representation of 



the SU(N) L x SU(N) R current (Kac-Moody) algebra (see for instance |TT[| ). 



6 Confinement 



In order to study confinement, we evaluate the Fredenhagen Marcu order pa- 
rameter ||12|| . It is defined by 



lim p(L) := lim E 

L^oo L— >oo ~T^, 
a, p=l 



N {L)\ N (L) 

(WW) 



(70) 



14 



where 

:= (ip^i-L, 0)U(C (L) )^\L, 0)) for arbitrary flavor a = 1, 2, ...TV , 

(71) 

and 

[/(£(£)) = e ie /c(i) ^l^ 1 " . (72) 
For the contour see Fig. la. denotes the Wilson loop 

W^ = e ieI W Mx)dX \ (73) 
along the contour Cffi given in Fig. lb. 



-L,0) 



(L,0) 



Fig. la Contour for the 
parallel transporter. 



2L 




Fig. lb Contour for the 
Wilson loop. 



The physical interpretation is simple: if there is no confinement and free 
'quarks' can be isolated, the free quark lives in a sector orthogonal to the 
vacuum sector (obtained by applying gauge invariant operators to the vac- 
uum); all 'mesonic' states belong to the vacuum sector. In particular p = 
follows in this case. In the case of confinement, however, one expects 'frag- 
mentation'. If one tries to separate a quark-antiquark pair, one is producing 
another quark-antiquark pair from the vacuum which binds to the separated 
pair and forms two mesons, i.e. a state in the vacuum sector. In this case one 
expects therefore p ^ 0. 

The computational trick is to rewrite the contour integral over as a 
scalar product with a current having its support on the contour, e.g. 

/ A,(x)dx, := f A,(x)jf L \x)d 2 x , (74) 

JC w J 

with 

n w(L) ( , _ ( 0( Xl + L)9(L - Xl ) [6(x 2 — L) — 5(x 2 + L)] \ ( , 
3 [X) ~ { 6{x 2 + L)6{L - x 2 ) [Six, + L) - 5( Xl - L)\ ) ' 1 <b) 
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(strictly speaking a limiting procedure is required here, which we leave as an 
exercise). Performing the functional integration one ends up with 

= J dficWe^^ = e -^ (L) > c ^ (L) ) . (76) 
For N a/ 3 in addition a propagator appears in the gauge field integral 

N$ = J d^[A] e ^^ (i) )G^((-L,0), (L,0);A) . (77) 

It is useful to rewrite the gauge field integration in terms of the field tp in- 
troduced in (16)-(18). In the chosen representation of the 7 matrices we have 

Nfff = and = n[^. Inspecting the quadratic forms that remains after 

the functional integration over ip one finds = N$ G JR. Putting things 
together gives 

= — — e l(L) (78) 

9 (2tt) 2 (2L) 2 ' 1 ] 

where the integral 1^ := 1^ + 1^ contains two parts 1^ and I^ L \ 1^ will 
be solved explicitely, and 1^ can be shown to go to a nonvanishing constant 
for L — > 00. I^ L > is given by (m 2 = e 2 N/n) 



To solve it, we transform it to polar coordinates and integrate over the angle. 
This leads to an integral over a Bessel function, which can (after some partial 
integration) be found in integral tables. The result is 



jvMWt +£• (80) 



where 7 denotes Euler's constant. The modified Bessel function K vanishes 
for L — > 00. The second integral 1^ is given by 

7r J p z + vn A y (Pi) 
4sin 2 ( Pl L)^[sin 2 (p 2 L) - sin 2 (p 2 L/2)]} = 

¥ / ^ V + ™ 2 )(P2) 2 COs(Pl2L) ( Sin2(p2L) ~ sin2 ^ 2L )} • ( 81 ) 
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Performing the p\ integration in the second term, one sees that this integral 
is exponentially decaying for L — > oo. The first term we integrate over pi first 
and find the cancellation of two terms linear in L. The remaining part gives 
an exponentially suppressed term, and the constant 1/N. 
Putting things together, we obtain for p(L) 



N 



(In the exponent we dropped terms that vanish for L — > oo.) For the confine- 
ment parameter this gives two different results depending on the number of 
flavors 

f J-^Ip 2 7+i for N — 1 
p = lim p(L) = 4 , e Ior iV - 1 . ( 83 ) 

L^oo rv ( for N > 1 v ; 

This result may appear surprising. Confinement of a quark-antiquark system 
depends on the number of flavors! (Up to a factor the N — 1 result can be 
found in 0.) 

Having solved the above case one can easily generalize the formula to prod- 
ucts of n spinors with pairwise distinct flavor quantum numbers. Define for 

n < N 

:= (^-£,0)...^(-L,0)f^^ . 

(84) 

The corresponding p^ n \L) is given by 

v /v( L ) f /v (L) 

P (L) - - (85) 

2 e ^ (L) +0(1/L) = - A_ (^^ e £(WJ_f ^ +0(1/L) 

(86) 

Performing the limit L — > oo we obtain 



p <») = l im p <») (i) = | f "* » = " . (87) 

for n < N 



Thus in the iV-flavor model an arrangement of iV 'quarks' is bound by a con- 
fining force to an arrangement of N 'antiquarks'. 

A further generalization shows that an operator of N + n, n < N quarks 
behaves like the product of only n quarks. The corresponding N is defined by 

N (L) 

^ ai ... aNa ' 1 ... a ' n (3 1 ...(3 N (3' 1 ...(3' n ■- 
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For the fields with flavors l,2,...n extra contractions are possible. Some of 
these cancel each other. The remaining terms have no dependence on (p and 
contribute only to the 1/L behaviour of the free propagator. The flavors 
n + 1, n + 2, ...N are seen to contribute a factor 

e (N-n ? m ? (g9) 

after performing the path integral. Collecting all terms one obtains 

p (N +n) = Um { N +n){L) = { const ^ foi :n = 0, N 

1 l^oo 1 I for < n < N v ' 

The physical interpretation suggested by this behavior of the Fredenhagen- 
Marcu order parameters is the following: The model has N distinct superse- 
lection sectors labeled by a charge Q that is defined only modulo N. To obtain 
a state in the the sector of charge Q = n, n < N, one applies an operator con- 
sisting of n 'quarks' and n antiquaries, separated by distance L and takes the 
limit L — > oo. 



7 Decomposition into clustering states and the 
vacuum angle 

The state (expectation functional) which we have constructed so far violates 
clustering, as we will show in this section. We will also construct an integral 
decomposition of this state into 'pure phases' labeled by a vacuum angle 9 and 
satisfying the cluster decomposition property. Thus we recover the structure 
that is usually associated with the existence of topological sectors of the gauge 
field. The words 'pure phase' should not be taken literally, because the different 
states are not related by a symmetry operation as in the case of spontaneous 
symmetry breaking. The transformations that intertwine the different ^-states 
- the axial U(l) transformations - are not symmetry transformations because 
the corresponding current is not conserved (due to the presence of the Adler- 
Bardeen anomaly), and therefore the corresponding Ward identities contain 
an anomaly term. 

7.1 Identification of operators that violate clustering 

To identify the operators that violate clustering, we start with an ansatz con- 
taining only the chiral densities ip^ P±ijj^ (P± :— (1 ±75) /2) and discuss the 
effect of adding vector currents later. Define 

N n a m a . . 

a=l i=l i=l 
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n ^ (a) {x't ] )P^ {a) (4 a) ) n ^ {a) (vf V-^ (a) (yt ] ) 



1=1 1=1 



n n^ (a) (^ (a) )^ (a) (4 a) ) n^ (a) fe (a V-^ (a) y a) ) 

a=l i=l i=l 



JV n' a m' a 

n n^ (a) « (a) )^ (a) « (a) ) n^ (a) (yf V-^ (a) (yi a) )) , (m) 

a=l i=l i=l 

where f is the vector of length r in 2-direction. Violation of the cluster property 
now manifests itself in a nonvanishing limit 

lim C(t) =:C^0. (92) 

T— >0O 

Because G a ^ a = 0, the second term contributes only if m a = n a ,m' a = n' a for 
all a — 1, ...N, and the first one only if 

n a + n' a = m a + m' a , a = l,...N. (93) 

After some reordering the first term reads 

N n a "i' a n' a m a 

n n^ a) (4 a) + r) n^ a) (yt } ) n#(^ a) ) n^ a) (^ a) +r) 

a=l 2=1 2=1 2=1 2=1 

m a TT-a m. a TIq 

nv4 a vs a) ) n^ 0) (^ o) + *) n4 a) (yi a) + *) n$Vi a) )> . (94) 

j=l i=l i=l j=l 

We dropped an overall sign depending on n a ,m a ,n' a ,m' a which is not relevant 
yet. Using the explicit form (20) of the propagator, one finds that the gauge 
field integral can be factorized off: 

(•) = /(-)o, (95) 

where (•) stands for the expectation value (94), and (-) denotes the same 
expectation value for e = 0. The factor I is given by 

i = Jd^ d w 2e ^ [s^-he:^ 

(96) 

It has the general structure 
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where we performed the functional integral and used (19) for the covariance 
C, to obtain after some reordering of the arguments 



V(x) = 4e 2 J 



d 2 p 1 — cos (pre) 
(2tx) 2 p 2 (p 2 + m 2 ) ' 



(98) 



Up to a trivial factor this is the integral defined in equation (79), which we 
already solved in the last section. The result is 



V(x) = ±Hx 2 ) + V'(x), V'(x) = l 



K I e\j—\x\ 



In 



7T 



2 V 7T 



+ 7 



(99) 

where we split V(x) into the logarithmic part and a part V'(x) which ap- 
proaches the constant j|[l n (f \Zv) +7] for x ^ oo. 

When applying (97) to (96), the sets {wi} and {zi} are given by 



€ {y\ b) + f, yT | I = 1, ...ma; k = 1, ...m' a , b = 1, ...N} , 



,/(<*) 



z f) e {x f) + fj x >w | / = l, ... na . k = 1, b = 1, ...TV} . 



/(«) 



(100) 



Inserting these arguments in (97) and collecting the terms with nontrivial r 
dependence, we obtain for the large-r behaviour of the gauge field integral a 
factor 



(Ea»»)(Ea "»-)+( E„"»«)(Ea"'0-(Ea»-)(Ea"-)-(Ea»»)(E.<) 



(101) 



The (-)o factor is 

AT . n a 



= n ( n (^ (a) + f ) n ^ (</r j ) n ^ « w ) n ^ or + *) 



a=l i=l 



V{y't ) )Xl^\x' { ? ) )W< {a)<Ja) 

i=l i=l i=l 



Tn a n a m a n a 

1=1 i=i i=i 1=1 

It can be expressed in terms of determinants. The general structure is 

n n 

c := ( n ^i(zi)ip 2 (wi) n v , i(^)^ 2 (^)) = 

i=i 

1 ■ 



i=i 



n 2 V . 1 

,- sign-—- - 



2tt/ 



XI si g n — 



1 



1 



„.( n) ^1 - ^tt(I) Z n — Wv(n) „./(„) ^1 - «V(1) - ^Tr'(n) 



(^) det (-4 



(103) 
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Determinants of this type can be rewritten using Cauchy's identity (see e.g. 
0) 

( 1 N ( ^!^ Ul<i<j<n{Zi ~ Zj){Wj ~ Wj) 

— n;: ; ..u,- -«-,•) — • (104) 

Hence we obtain 

n Hl<i<j<n( z i ~ z j) 2 ( w i — w j) 2 / 1nt -\ 

Ca - " rewfe^) 5 ' (105) 

Inserting the sets (100) for w^Zi and collecting the r dependent powers gives 
the large r behaviour 

(t J = [ T ) ■ ( 106 ) 

0=1 

Combining this with (101) gives the total large-r dependence for the first term 
of C(t) 

i /„ „ \r„i ™'\ Y^ 



' 1 \ i E„,i P i( ,1 »- m «)( n r m i)-E^i(' l «- m «)(''«- m Ii) / 1 \ E 



Violation of clustering can only appear if the exponent vanishes and C(r) 
cannot approach zero. Inserting the condition (93) we can rewrite the exponent 

TV 1 N 

E = J2( n a-m a )(n a -m a ) - — ^ {n a -m a )(n b -m b ) = 

0=1 ^ V a,6=l 

1 * 

T7 {n a -m a )R ab (n b -m b ) , (108) 

iV o,6=l 

where the matrix R was already defined in (25), section 3. There we solved 
the corresponding eigenvalue problem. We found one eigenvalue 0, and iV — 
1 eigenvalues N. The eigenvector x° to the eigenvalue is given by x° = 
l/v~N(l, 1, ...1) T . Hence the quadratic form x T Rx is positive semidefinite, 
and vanishes only if x is a multiple of x°. This implies that the exponent E is 
nonnegative and vanishes only for 

n a — m a = m ' a ~ n 'a = U i ™ E 7L . (109) 

In the special case n = the second term in (91) does not vanish. It cancels 
the contribution of the first term. Hence we find no violation of clustering 
for n = 0. As already mentioned, for n ^ the second term vanishes due to 
G aa = 0, and no such cancellation is possible. 

How does this picture change when we allow vector currents as well? First 
we notice that vector currents do not contribute to the gauge field integral. 
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Consider e.g. the term ^^(x)^' (x). The ^'(x) generates a propagator 
G 2 i(-, x), the ip^ix) a propagator G 2 i(x, ■). Inspecting (20) immediately shows 
the cancellation of the ip dependence in the product of the propagators. Hence 
each vector current can only contribute a 1/r from the free propagator. 

Nonvanishing results remain only if the flavors that occur in the vector 
currents can contract entirely. Thus we have to consider only 'closed cycles' 
like e.g. 

(110) 



,(2) 



t(1), 



^ {1 \x)^(x) ^\y) lv ^\y) ^\z) lu} ^\z) . 

In principle there are two possibilities to distribute the space-time arguments 
x, y, z. If they are all in one cluster they do not bring in any r dependence. 
They do not modify the clustering, only the constant C. If one distributes the 
closed cycle over both clusters then the situation changes. Each vector current 
with a partner in the other cluster contributes a 1/r from the free propagator. 
This implies that any combination of vector currents alone clusters. Never- 
theless a combination of vector currents together with the ansatz (91) could 
violate clustering. But we saw that the gauge field integral contributions from 
the chiral charges %l)P±%p can at most compensate the 1/r of these charges, 
nothing else (x T Rx is positive semi definite, compare (107)). Hence adding 
vector currents that can only contract between the clusters can at most en- 
force operators to cluster, never create extra powers of r that lead to violation 
of clustering. The same is true when inserting currents of the Cartan subal- 
gebra where one has to introduce a point splitting regulator. The gauge field 
transporter that connects ip(x — e) and ip(x + e) only is a modification within 
a cluster that does not change the clustering behaviour. 

We now want to analyze the symmetry properties of the nonclustering 
operators. To make our notation more convenient we introduce 



N 



o±({x}) -.= n^ (a) (^ {a) )^±^ (a) (^ (a) ) 



(111) 



It will turn out that the lack of clustering of 0± is related to the fact that 
they are singlets under the conserved symmetry group U{l)y x SU(N) L x 
SU(N)r, but transform nontrivially under the explicitly broken U(1)a- To 
obtain operators that transform under a single irreducible representation of 
the symmetry group, we antisymmetrize 0± with respect to the flavor indices 
and call the result 0±. 



N(N-l) 

Ol({x}) : = (-l)-HH 



N 



a I 2 



N 



^Esign^n^V" ) 

tt' a'=l 2 



(112) 
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The global sign comes from shifting all ip to the left. Using 

N N 

n = ^(vr- 1 ) n ?r (a V o) ) , (us) 

a=l a=l 

one can express 0± as well in terms of symmetrized space time arguments, 

1 N _^ 

±(M) : = 7W E II ^ a \x^)P±^ a \x^) . (114) 

^ iV 'J 7r,7r'a=l 

Since the constant C := lim r _ >00 C(r) is invariant under the permutation of 
arguments within a cluster (compare (97), (104)), the latter expression shows 
explicitely that the constant C is the same for 0± and 0±. From (112) one 
easily reads off the invariance of 0± under 

U(l) v x SU(N) L x SU(N) R . (115) 

We end up with the following picture of the clustering problem: The prototype 
of a correlation function that violates clustering is given by 

tl y~n n/ to! 

C(r) = ( II O a + ({*i + r}) II W-dVi + r}) U O a + ({x[}) J] Ot{{^))- 

i=l i=l i=l i=l 

Tl Tfl Tl! m/ 

n^(M)n^(M)> <n^(K})n^({^})>> 

i=l i=l i=l i=l 

with the condition 

n-m = -ri + rri E ^\{0} . (117) 

Insertion of closed cycles of vector currents into a cluster does not change 
the clustering behaviour. If one inserts vector currents that can contract only 
to a partner in the other cluster, operators that violated clustering before 
now become operators obeying the cluster property. Of course it is possible 
to generalize the operators 0± further by e.g. splitting the arguments and 
connect them with a parallel transporter. Since this is a modification within 
a cluster, the extra terms in the functional integral will not depend on r and 
only modify the constant. 

For completeness we quote the constant C = lim T ^ 00 C(r) with C{r) de- 
fined in (116) 

C = F({ Xi },{ yi }) FiWMlti)- (118) 
The dependence on the space time arguments factorizes into two parts that 
depend on the arguments in the two clusters. This function T is unique for an 
operator, if only the other operator has the right quantum numbers to form a 
pair that violates clustering. It is given by 

lll<i<j<n \ x i x j J lll<i<j<m !Jj J 



N 



r({zih{vi}) = U 



a=i 



n?=i ut=i (4 a) - yf) 
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E^(* ( (a) -»f ) )-3i;v(xW-*5 6) )-jEv(»J- ) -»f) ^j_y 



N(n+m) 



e 2 N 
An 



,2-> 



(119) 

E denotes summation over all possible values of the indices a,b,i,j. This 
factorization property remains valid when vector currents are inserted. 



7.2 Decomposition into clustering states 

To decompose the vacuum state of the theory in terms of clustering 9 vacua 
we use the charge that is associated to the axial U(l) transformation 

An arbitrary product B of tfj^\x), ip^, \x') transforms under U(1)a 

B{{x}) — > e tm B({x}) , meTL. (121) 

(More generally we can consider observables that are sums of operators with 
definite transformation properties under U{1)a)- Define the corresponding 
charge Qs(b) as m. Obviously 0± have the charge ±2N. We will now de- 
compose the expectation functional (•) into states (-)g labeled by a parameter 
9 E [— 7r, 7r] defined as follows 

(B({*})) 9 ■= e l6Q ^ }}m(U T (B) B( W )> . (122) 

The set of 'test operators' U t (b) is defined by 

U (b) = { ^ (n) (W)nr=i^({^}) for Qb(b) = ±2niV,n > 1 , 
T ^ ' | 1 otherwise . 

Up to the requirement of being nondegenerate, the arguments {x} are arbi- 
trary. The normalizing factor J\f( n \{x}) is defined such that 
lim T /^ 00 (W T /(et) U t (b)) = 1. It can be read off from (119) 

n li (4 a M ) 1 3 > . 

a=l l<i<j<n 

(124) 

The expectation functional (states) (-)g have the following properties: 
Theorem 3 : 

i) The state (•) constructed initially is recovered by averaging over 9 

( • ) = ^ J J ■ )ed9 . (125) 

ii) The cluster decomposition property holds. 



/ l y Nn 




V2W 


47T 
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Proof: 

i) : The averaging procedure leaves a nonvanishing result only for operators B 
with vanishing charge Q*,(b) = 

(B) = if (B)edd = -L f e ^»W dfl lim (W t (b) B> = 5q 5(s) , (B> . 

Z7T J -IT Z7T T^OO 

(126) 

In the last step we used U t (b) = 1 if Qb(b) = 0. To complete the proof, one 
has to check that in the state (•) constructed initially the vacuum expectation 
values of operators B with Q$(b) ^ vanish (cf. also In our representa- 

tion of the 7 matrices this can be seen immediately. Q${b) ^ means that the 
number of ipi is not equal to the number if) 2 i i ) 2- Since G aa = 0, for each ip\ 
there has to be a if) 2 an d f° r each ip2 a ip 1 to give a nonvanishing contribution. 
But this implies that the number of fields ipi^i m & is equal to the number 
of fields ip 2 i4'2- Hence 

(B) = for Q 5 {b) ^ . (127) 

ii) : Let A and B be arbitrary operators. Define 

C e {r) := hm [(A(r)B(0)) e - (A(0))e (B(0)) e ) . (128) 

We do not display the dependence on the space time arguments {x} explicitely, 
only the dependence on the shift variable r. Depending on the axial charge 
Qs of the operators A, B, AB, we have to insert the different definitions of 
(■)g. We introduce the following convenient notation: An operator for which 
the first alternative in equation (123) holds we call a 'type I' operator, the 
operators where the second alternative holds is called 'type II'. We remark 
that all operators of type II cannot have the problem of violating the clus- 
tering condition. Even among the type I operators there are examples that 
are not able to violate clustering (like e.g. (ip ° P + ip^) N which does not con- 
tain a SU(N)l x SU(N)r singlet part). The operators with the structure 
IT 0+ IT O- EI JP , n-m^0 we call 'type V for violating. 
We have to distinguish the following cases 



case # 




Qb{b) 


Q 5 {ab) 


1 


II 


II 


II 


2 


II #0) 


II #0) 


I 


3 


II #0) 


I 


II 


4 


II (=0) 


I 


I 


5 


i (=<?) 


I (="<?) 


II (=0) 


6 


i 


I 


I 



In brackets (•) we denoted facts that necessarily follow. For example in the 
case 2 the requirement Q§{a) G II, Q 5 (b) G II, Q$(ab) G I implies Qs(a) ^ 
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and Qs(b) ^ 0. If one of these two had charge zero, the other operator would 
be of type I, since Q$(ab) = Q<s(a) + Qs{b). 

Case 1: 

C e {r) = (A(r)B) - (A) (B) ™ , (129) 

since if A and B are of type II, they do not form a pair that violates clustering. 
Case 2: 

C (t)= lim (U t ,(ab)A(t)B) - {A){B) T ^ 

t'— +oo 

(A) lim (U t ,{ab)B) - (A)(B) = 0, (130) 

t'— >oo 

where we used the fact that (A) factorizes, since A is type II and (A) = since 
Qs(a) 7^ (compare the note in the table). The interchange of the r, r' limits 
is justified, since all the functions involved are continuous in these variables 
and bounded (the exponent E in equation (107) is nonnegative). 
Case 3: 

C e {r) = (A(t)B) - (A) lim (U t ,{b)B) ™ , (131) 

for the same reasons as in the last case. 
Case 4: 

C e {r) = lim (U t ,(ab)A(t)B) - (A) lim (U t ,(b)B) ™ 

t'— >oo t'— >oo 

(A) lim (U t ,(ab)B) - (A) lim (U t >(b)B) = . (132) 

t'^oo t'— »oo 

Again (A) factorizes, and U t >(ab) = U t >(b) since Q 5 (ab) = Q$(b). 
Case 5: 

C q (t) = (A(t)B) - lim {Ut>(a)A) lim (U t „(b)B) ™ 



Fa^b — FaFb — for .4., £> of V-type , 
otherwise . 



(133) 



In the first case we used the factorization of the argument function jF^g intro- 
duced in (118). 
Case 6: 

C e (r)= lim (U t ,{ab)A{t)B) - lim (U t/ (a)A) lim (Ut»(b)B) ™ 

t'^oo t'^oo t"^oo 

f ^4 linv-oo ^/ t ,b - ^U^s = for A, B of V-type , 

1 otherwise . ^ ' 

To justify the first case, one has to show that linv^oo Tu t ,b — ?b- This can 
be seen immediately from (119). Shifting r' in Tu ,b corresponds to shifting 
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the arguments of the test operator U t i{b). For example this could be the set 
(refering to (119) for J 7 ) 



f (a) I / 

l n 



<i<n;a = 1,2,..JV} , (135) 



where n' < n depends on the charge Q^{b). The t' terms cancel for r' — > oo, 
and what remains is J-jg, since the normalization of U t >(b) cancels exactly the 
extra terms. □ 

This concludes our analysis of the vacuum structure of the iV flavor Schwinger 
model. The result is in full agreement with the picture that is conventionally 
deduced from the discussions of topological sectors [|l(J , ]17 . 
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